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VARIANCE DETECTOR REPORT 



INTRODUCTION 

This report documents the program "vardet.c" written in July, 1984 for 
post- (FFT) processor activity/gap detection. The algorithm used is based on 
the assumption that the variance of the logs of the frequency- domain magni- 
tudes of a signal is greater than the corresponding variance of the channel 
noise. 

The noise is assumed to be quasi-ergodic, that is ensemble and single FFT 
averages and variances are assumed to be the same, except that these measures 

are allowed to vary slowly in time using a weighted, recursive technique. 

Trade-offs between Type I (false detection) and Type II (missed signal) 
errors can be controlled with an interactively defined program parameter. 

The program is designed to warm-up gracefully, if an a priori estimate of an 

upper bound on the variance of the logs of the noise magnitudes is available. 

In the absence of such an estimate, the warm-up period can be extended inter- 
actively to allow the program to learn the system noise levels. 

Additionally, this program allows the addition of noise to an input signal 
so that program parameters may be tested under various signal-to-noise ratios. 

The output of the program is in two forms. Indication of activity or gap is 
printed to a file called "GRAF" for graphing on the line printer. For periods 
in which activity is detected, the raw, pre-FFT data is printed to a file 
called "ACTIVE" for D/A conversion to recover the gap-surpressed signal. An 
estimation of the signal-to-noise ratio is also available for printing. 



TASK DEFINITION 



The problem is to design an activity/gap detector for digitized signal data 
of the following three types: speech, fsk, and psk. Running time is limited 

by the specifications of the mode analysis unit into which the detector is to 
be incorporated, and number of operations is limited by hardware considerations. 
The measure of effectiveness can be given by the complement of the following 
function : 



1) t = C * 




MISSED DETECTIONS 



d(snr) 




FALSE ALARMS 



d(snr ) 



where snr is the signal-to-noise ratio and the limits of integration are over 
the expected operating range, probably 5 to °° dB. This equation implies that 
missed detections (Type II errors) and false alarms (Type I errors) will vary 
with snr, which is expected. The constant C implies that missed detections 
and false alarms might be weighted differently when measuring effectiveness. 

We have assumed that missed detections are considerably more serious than false 
al arms. 



The task then, in summary, is to develop an alogrithm for signal/gap detection 
that minimizes 0, which meeting the constraints of running time and implementing 
hardware. 



THE APPROACH 



It was noticed upon examination of typical speech and fsk data that the variance 
of the logrithms of the FFT magnitudes increases when signals are present in a 
channel. A detector designed to exploit this finding would have several desir- 
able characteristics. First, variance is easy to compute, requiring only the 
summing then squaring and squaring then summing of values (about 3n operations). 
Second, the variance seems to be a fairly sensitive indicator of activity, much 
more sensitive than an energy level (mean value) detector. Third, use of the 
logrithm of the magnitude of the FFT makes the algorithm independent of time- 
invariant system gain levels. Note that: 

log(ax) = log a + log x = constant + log x 

But when considering the variance, 

variance [constant + f(x)] = variance [f(x)] 

This implies that 

2) variance Clog(ax)] = variance [constant + log x] = variance [log x] 

Thus demonstrating the algorithm's independence of any constant gain a. 

The use of logs complicates the algorithm, increasing its running time and pre- 
venting the calculation of the measure variance [log || FFT(signal + noise) ||] 
as a function of the signal and noise distribution parameters. 



THE ALGORITHM 

The variance-based detection algorithm can be given as: 

If variance [log I I FFT(waveform) | | ] > some threshold 
declare activity. 

Otherwise, declare inactivity. 

The problem becomes establishing an appropriate threshold. Clearly, a low 
threshold will lead to a high false alarm rate, while a high threshold will 
cause missed detections. The threshold must be related to the variance of 
the noise alone as follows: 

threshold > variance [log I I FFT(noise) | I ] 

Hereafter, the right hand side above will be referred to as var(noise) for 
simpl ici ty. 

The communication channels under consideration are 4 kHz wide, but only the 
300 - 3700 Hz band contains activity. This means that the first and last 
300 Hz bands in the channel should contain only noise. So with a 64 point 
FFT, the middle 56 bins will contain the activity, if present and the first 
and last 4 bins will contain the FFT of the noise only. In practice, we have 
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noticed that the two bins immediately adjacent to the activity band on each 
side are contaminated due to filter and windowing imprecision. This leaves a 
total of 4 bins, 2 above and 2 below the activity channel, as measures of the 
noise only. 

Four samples, however, are too few to get an adequate measure of the noise 
variance. If we can assume the noise to be ergodic, we can use the variance 
over time as a measure of the variance over frequency. That is, if the noise 
is ergodic, the statistical measures at any time over all frequencies will be 
the same as the measures at any frequency over all time. Because we have 4 
noise bins, we will have 4 independent measures of noise variance over time. 

Given the assumption of ergodicity, these will be close to the variance of the 
noise over frequency and can be used in setting the threshold. 

We experimented with different ways of using these 4 measures to arrive at a 
single estimate of var(noise). One method used, but discarded, is to take the 
maximum of the 4 values at any given time. The problem with this method is 
the extreme variability between successive threshold values. A second method, 
found to be more stable, is to take the average of these 4 variances as an 
estimate of the true variance. To establish a threshold, we increase this 
value by multiplying it by a "kluge" factor (K) to be determined empirically 
and depending upon the distribution of the noise over time and the trade-off 
between Type I and Type II errors. Thus, 

3) threshold = K * average over 4 bins Cvar(noise) over time] > var(noise) 

We can summarize the detection algorithm as: 

If variance across inner 56 activity bins at any time [log | | FFT(waveform) | | ] 

is greater than 

K * average across the 4 outer bins {variance over time of each of the 4 bins 
[log | IFFT(waveform) | |] } 



declare activity. 

Note again that the variance across the inner activity bins is at a single time 
window, while the variance of the noise is measured across all time windows. 

In application, we are using a 64 point FFT requiring 128 data points. In 
general, activity will not truely begin or end at the start of one of these FFT 
cycles. Therefore, we expect a one cycle uncertainty regarding the true start 
and end of activity. A data sampling rate of 8000 for the 4000 Hz transmission 
channel requires 62.5 FFT cycles per second of 128 data points for real time 
processing. Thus, one cycle would translate to a 16 millisecond uncertainty 
regarding the true start and end of activity. Speech activity often begins and 
ends at low energy and, presumably, low variance levels, so the true uncertainty 
in practice would be even higher. 



ALLOWING FOR TIME-VARYING NOISE DISTRIBUTIONS 

We would like to make as few assumptions regarding the noise distribution as 
possible. We have thus far assumed only ergodicity. It may be possible to 
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create an algorithm which is robust against violations of this assumption by 
allowing for a slowly (with respect to the time of one FFT) time-varying noise 
distribution. This could be accomplished by using a moving, weighted measure 
of the variance across time of each of the 4 noise bins. 

We will start by finding the moving, weighted mean of the magnitude of the FFT 
for each of the 4 noise bins. 



a°x + a x x . + a 2 x + a n x 0 

x _ n n-i n -2 

n a 0 + a 1 + a 2 a" 

where x is the magnitude of the FFT of the i th bin at each time window and 
o < a < 1 is the attenuation constant. 

For the purposes of calculation, this can be written recursively as 



5) 



X 

n+1 



a. 

a*X 




+ 



x 

n 



d = a*d + 1 

n+1 n 



Now we calculate 
6) VAR = 



the moving, weighted variance based on this weighted mean. 

a ° (x n " V 2 + al(x n-i ' X n-i )2 + a " (x ° " X ° )2 

a 0 + a 1 + a 2 + a" 



This can be written recursively as 



7) 


U 


a*U + (x - X ) 2 




n+1 


n n+1 n+1 




d 


a*d + 1 




n+1 


n 






U 




VAR 


= n+1 




n+1 


(n+1 )*d 






n+1 


The 


use of a 


moving, weighted variance 



however, has effects not anticipated. 
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Basing the variance on a moving mean biases the measure, causing it to be an 
underestimate. There is a mechanical analogy illustrating this effect. If a 
pendulum is allowed to swing on a moving pivot, the pendulum will swing faster. 
The pivot moves in the direction of the center of gravity (moving mean). The 
faster swing indicates that the momemt of angular momentum (second moment or 
variance) has been lowered. 

Even if the noise distribution were completely specified, finding the degree 
to which the variance is underestimated by this method would still be a for- 
midable task. Recall from equation 3), however, that the average of the 4 bin 
variances will be multiplied by the kluge factor K. If the noise distribution 
doesn't vary too wildly, the kluge factor can be turned empirically to com- 
pensate for the underestimate of the variance. 

It will be noted that the denominator in both the moving mean and variance cal- 
culations is a geometric series with a constant 0 < a < 1. This series converges 
to l/(l-a). For the purposes of implementing this algorithm into hardware, this 
limiting value is used through all calculations. This, of course, means that 
for the early cycles the mean will be grossly underestimated and the variance 
will be grossly overestimated. 



GRACEFUL WARMUP 

It was mentioned above that approximating the denominator of the moving statis- 
tics by the limit of the geometric series causes early measures to be very poor. 
There is a second reason to expect problems in the early cycles. 

Note that in equation 5) the last term on the right hand side (x 0 - X 0 ) is 
always zero because 



0 1 

This means that VAR = 0 and, by equation 3), the threshold = 0. For small 
sample sizes (shortly after the algorithm is turned on) the variance, and hence 
the threshold, will be grossly underestimated. Consequently, early waveforms 
will invariably declared activity, whether or not activity is truely present. 

This is compensated for in two ways. First, the algorithm can be given an 
adequate chance to warm up. In the computer program to be presented in following 
sections, the warm up period is interactively specified by the user. Second, 
if an upper bound on var (noise) is available a priori, it can be used to replace 
(x - X) in equation 6). The following program also incorporates this method. 



NOISE TRANSIENT AND ACTIVITY DROP-OUT FILTERING 

We noticed in application that in the midst of a string of non-activity de- 
clarations, a spurious declaration of activity would often take place. In 
actuality, speech activity of single cycle duration would be most unusual. 
Similarly, single cycle drop-outs would often occur in the midst of a string 
of activity. To eliminate these problems, we adopted a "2 out of 3" filtering 
criterion. The filter is activated at the discretion of the program user by 
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interactively responding to query regarding filtering. The filter simply allows 
activity to be declared only if 2 of the last 3 cycles indicated activity. 
Similarly, non-activity is declared only if 2 of the last 3 cycles indicated 
non-activity. 

One effect of the filter is to cause the first cycle of a string of activity to 
be invariably declared as non-activity. If real-time delay and storage require- 
ments are not a problem, all data may be saved through the following cycle. If 
a cycle is the first in a string, the data may be retrieved after the next cycle 
and then written to the output files. 



ESTABLISHING OPTIMAL PARAMETER VALUES 

Parameters K, a, and the initial value of the var(noise) can be optimally set 
by exhaustion using equation 1). Due to the magnitude of this project, it has 
not yet been done. Rather, preliminary values of these parameters have been 
found which yield satisfactory results for the waveforms thus far tested. 

Under the previous assumptions, the value of K should be greater than 1. 
Experimental evidence, however, has lead us to believe that for signal and noise 
environments thus far tested, K should be less than 1 for optimal performance. 
There are two reasons for this surprising result. The first is the judgment 
that Type II errors (missed detections) are more serious than Type I errors 
(false alarms). Recall that lowering the threshold decreases missed detections 
at the cost of false alarms. The second reason, however, is fundamentally 
related to the characteristics of the noise. We noticed that when activity is 
present, the magnitude of FFT's of the 4 outer bins drops about 40% for the 
waveforms tested. This would be true if an automatic gain control were present 
somewhere in the system. A time-varying gain causes the noise to be non-ergodic. 
The variances over time of the outer bin FFT's increase because the magnitudes 
change (decrease) whenever activity is present. All this serves to overestimate 
the noise variance. Thus, it appears that an optimal value of K will be some- 
what below 1. In the following program, the user specifies K interactively. 

The attenuation constant must also be specified, by nature, between 0 and 1. 

A low value allows greater movement of the weighted mean, thus lowering the 
measure of the variance and causing a high false alarm rate, which can be com- 
pensated for by raising the value of K. A high value causes the mean and variance 
to become more insensitive to changes in the noise parameters. This may be 
desirable as it causes the threshold to remain stable through brief periods of 
noise instability. 

The parameters were chosen on the basis of tests on a 102400 value data set. 

This set yielded 800 FFT cycles. Approximately 258 of these cycles were 
identified ahead of time as containing activity. The remaining 552 were felt 
to be non-active. We found that for the single data set tested, with no added 
noise, the values K = 0.8, a = 0.997 and U = 57.0 worked well, yielding roughly 
41 missed detections and 26 false alarms. For the purposes of speech compres- 
sion, this result is entirely satisfactory. 



THE DETECTION PROGRAM 

A flow chart of the program is given as Figure 1. The program itself, written 
in the C programming language, is included in the appendix along with a variable 
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list. The main program requires four external functions (subroutines). These 
are: iread and iwrite, both written by Steven D. Beck to read into the main 

program the digitized signal values and to write to an output file the gap sur- 
passed digitized signal; fft2, provided by Steven Beck, to perform the Fast 
Fourier Transform; and actgraf2, written by Walter Underwood to create a file 
for graphing on the line printer the indication of gap or activity. 



NOISE SIMULATION 

The program has the capability of adding noise to the signal prior to the 
variance calculations so that program performance may be tested at low signal - 
to-noise ratios. The program assumes that the warmup period consists of non- 
activity. The complex FFT data for each warmup cycle is stored. We exploit the 
linearity of Fourier transforms to add the warmup FFTs multiplied by an inter- 
actively set gain to the FFTs of subsequent cycles to get the equivalent FFTs 
of noisy data. The warmup FFTs are added to the subsequent FFTs cyclically. 

The advantage of this method is that no assumptions regarding noise, except 
stationari ty, are made. We are assuming, however, that the warmup period 
contains noise only. Consequently, the noise simulation works only with data 
files that begin with non-activity. 



ALGORITHM PERFORMANCE 

As previously mentioned, parameter values were found which yielded approximately 
a 16% missed detection rate and a 10% false alarm rate, basing both percentages 
on the number of cycles actually felt to contain activity. These values were 
found with data exhibiting approximately a 20 dB signal -to-noise ratio for single 
activity cycles. Added noise caused graceful performance degradation, primarily 
increasing the false alarm rate. 

An obvious weakness of the algorithm, however, was its inability to detect the 
spoken word "six" as activity. The "s" and "x" sibilate sounds are broadband. 

The effect of these sounds is identical to adding noise to the activity band. 

The variance of the FFT of this band is not increased by the presence of these 
sounds, hence their undetectabil i ty. This appears to be a fundamental weakness 
of the algorithm and not the result of poorly chosen detection parameters. 



HARDWARE IMPLEMENTATION 

A flow chart indicating possible hardware implementation is attached. 



SUMMARY 

A detection algorithm has been written that gives good results for speech 
activity and also is simple enough for hardware implementation. Performance 
with fsk and psk signals has not yet been tested, but we expect satisfactory 
performance with these signal types as well. No assumptions regarding the 
distribution of the underlying noise are made, although alogrithm parameters 
are tuned on the basis of empirical knowledge of the noise. The system is ro- 
bust against changing noise parameters and degrades gracefully with decreasing 
signal-to-noise ratio. The algorithm has proven unable to detect sibilates, 
but otherwise gives excellent performance. 
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APPENDIX 



The following appendix contains the flow chart, the variable list and a listing 
of the program "vardet.c". 
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VARIABLE LIST 



ang 


added noise gain, used when simulating added noise 


a ten 


atenuation constant for moving weighted mean and 
variance. 0 < aten < 1 


avesnr 


average signal -to-noise ratio over 6 declarations of 
activity 


avevar 


average of var(noise) across the s4 outer bins 


counter 


index used when windowing digitized waveform prior 
to FFT 


countO 

countl 

count2 


flags used for 2 out of 3 filtering. l=activity, 
0=inactivity. CountO is the flag for the most recent 
FFT, countl for the previous, count2 for second 
previous . 


denom 


denominator in weighted, moving mean and variance. 
Geometric series approximated by limit value 1/(1 -a ) . 


df 


degrees of freedom of inner bin variance = # of 
bins - 1 


end 


flag transferred from iread(). 1= end of data 


fin 


general purpose character used for file names 


fP 

fph 

fpw 


file pointers 


h[l 28] 


coefficients for windowing waveform prior to FFT. 
read from file "FFTWIN" . 


i 


general counter 


idf 


1 /df 


im 


temporary variable name of imaginary part of waveform 
FFT plus imaginary part of added noise FFT 


inoise[j][64] 


imaginary part of the FFT of the noise, j th noise 
sample, i th frequency bin. IT IS ASSUMED THAT THE 
WARMUP PERIOD CONSISTS ENTIRELY OF NOISE. 


irec 


counter used exclusively by iwrite () 


isactive 


indicator transferred to actgraf2() for graphing. 
Inactivity, 0=inactivity 


k 


counter for number of FFTs performed 
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klug 


klug value used to multiply variance to establish 
threshol d 


mag[64] 


magnitude of the FFT of the waveform plus added noise 


n 


number of data points in the FFT (128) 


nar 


interactively set flag to indicate whether continuous 
display of program variables k, varib, avevar, thresh 
and snr is desired. l=desired. 


nisact 


counter for actual number of declarations of activity 


nn 


number of cycles indicating inactivity 


ns 


number of cycles indicating activity 


r[64] 


real part of the waveform FFT transferred from fft2() 


re 


temporary variable for real part of waveform FFT 
plus the real part of the added noise FFT 


rnoise[j][64] 


real part of the FFT of the noise, j th noise sample, 
i th frequency bin. IT IS ASSUMED THAT THE WARMUP 
PERIOD CONSISTS ENTIRELY OF NOISE. 


snr 


si gnal -to-noise ratio, calculated for each FFT 


ss 


sum of the squares of the logs of the magnitudes of the 
inner 56 activity bins, used to calculate the variance 


ssn 


energy in the outer 4 bins, used for calculating signal 
to-noise ratio 


sss 


energy in the inner 56 activity bins, used for cal- 
culating signal -to-noise ratio 


sum 


sum of the logs of the magnitudes of the inner 56 
activity bins, used to calculate the variance 


suma[4] 


the numerator in equation 5), used to calculate the 
moving, weighted mean 


temp 


log base 10 of the magnitude of the FFT bin currently 
under consideration. 


thresh 


threshold value for declaration of activity 


two 


interactively set flag to indicate whether 2 out of 
3 criterion for activity declaration is desired. 

1 =desi red. 


Unext 


U given in equation 7) used for calculating the 

n+1 

moving variance. 
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Uold U given in equation 7) used for calculating the 

n 

moving variance. 

varib variance of the inner 56 frequency bins 

warmup interactively set number of warmup FFTs 

wtave weighted, moving mean used to calculate moving 

variance 

x[128] digitized waveform data, read from input file 

y[64] imaginary part of the waveform FFT transferred from 

fft2() 
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y'*********************** PROGRAM: vardet . c ************************ 

* 

* Written July, 1984 by Dr. J.L Wayman to do activity/gap detection. 

* 

* Program is written in the C language. 

* 

* Subroutines Required : iread () , fft2 () , actgraf2 () , iwrite () 

* 

V 

#include <stdio.h> 

#include <math.h> 

ttdefine PMODE 0644 



short int fpr, fpw; 
main () 

{ 



float 


r [256] , 
y [256] , 
h[128]; 


double 


moise [16] [64] , inoise [16] [64] , 
mag [64] , suma [4] , 

Uold [4] , Unext [4] , 
sum, re, im, avevar , 


V 


ss, avesnr, 
sss, ssn, snr, 
idf , df , ang, aten, 
denom, wtave , klug, 
varib, 

temp, thresh; 


long int 


two, nisact , 
nn, ns; 


short int 


x[128] , 
n, 

counter, 
end, warmup, 
i, j , k, 

countO , countl , 
count 2, 
nar , 
ir ec ; 


int 


isactive; 


char 


fin [20] ; 


FILE *fopen(). 


* fph; 


FILE *fopen(). 


*fp; 
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^/************************* OPEN X/O FILES ****************************y/ 



if ( (fp = fopen ("GRAF" , "w") ) ==NULL) 

< 

printf ("Error cannot open GRAE\n") ; 
exit () ; 

> 

i f ( ( fph = fopen ("EFTWIN" ,. "r") ) ==NULL) 

printf ("\nERROR: cannot open %s for reading\n" , "EFTWIN") ; 
exit () ; 

). /* Read coefficients for 

for (counter=0; counter<128; counter++) /* EFT windowing 

{ 

fscanf (fph, "%f ", <Sh [counter] ) ; 

> 

printf ("\nEnter output file name:"); 
scanf ("%s", fin) ; 

if((fpw = creat (fin, PMODE) ) =-l) 

{ 

printf ("ERROR, cannot open output file for writing\n") ; 
exit () ; 

> 

printf ("\nEnter input filename : ") ; 
scanf ("%s", fin) ; 

if ( (fpr=open (fin, 0) ) = -1) 

printf ("\nERROR: unable to open %s for reading\n" , fin) ; 
exit () ,-i 
> 



y* ****************** ★ INTERACTIVELY SET PARAMETERS ********************/ 

printf ("\nEnter added noise gain:"); 
scanf ("% f " , &ang) ; 

printf ("\nEn ter attenuation factor for moving averages:"); 
scanf ("%f" , &aten) ; 

printf ("\nEnter kluge factor for threshold:"); 
scanf ("%f",&klug) ; 

printf ("\nEnter warmup period (<=16) :"); 
scanf ("%d" , Awarmup) ; 

printf ("\nNarrative mode (l=yes , 0=no) ?") ; 
scanf ("%d" , &nar) ; 

printf ("\nTwo out of three filtering (l=yes, 0=no) ?") ; 
scanf ("%d" , &two) ; 
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*********************** * INITIALIZE VARIABLES *********************/ 

n = 128; 

irec=0; 

df = 55.0; 

end = O; 

idf = l.O/df; 

k = O; 

ns=0; 

nn=0; 

nisact=0; 

denom =1/(1 -a ten) ; 

isactive=0;. 

countO=0; 

countl=0; 

count 2=0; 

for (i=0; i<4; i++) 

{ 

suma [i] =0; 

Uold [i] =57.0; 

> 



/******************** READ AND PROCESS WAVEFORM DATA *****************/ 

while (1) /* While there is still raw */ 

{ /* waveform data to be read */ 

iread (x„ n, &end) ; 

if (end = 1) break; /* read and window data */ 

for (counter=0; counter <n; counter++) 

{ 

r [counter] = x [counter] *h [counter] ; 
y [counter] = 0.0; 

}' 



fft2 (r,y, 128, 7, 1) ; 

i f (k<warmup) 

{ 

for (i=0; i<64; i++) 

moise[k] [i] = r[i]; 
inoise [k] [i] = y[i]; 

> 



/* Go to subroutine for EFT */ 



/* During warmup period */ 
/* store away noise for later */ 
/* use in simulation */ 



for (i=0; i<64; i++) 

{ 

j=k%warmup; 

re = r[i] + ang*rnoise [j] [i] ; /* Add noise to signal and */ 

im = y[i] + ang* inoise [j] [i] ; /* take magnitude */ 

mag[i] = sqrt ( (double) (re*re +im*im)); 

> 



sss = O; 
ssn = O; 



/* These variables will be used */ 
/* to establish snr for this cycle*/ 
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/************************ COMPUTE ENSEMBLE NOISE VARIANCE *****************/ 

for (i=0; i<2; i++) 

{ 

temp=loglO (mag[i] ) ; /* First two FFT bins */ 

suma[i] = aten*suma [i] + temp; 
wtave= suma [i] /denom; 

Unext[i] = Uold[i] + (temp*temp -2 .0*temp*wtave + wtave*wtave) ; 
ssn = ssn + mag[i] *mag[i] ; 

temp = loglO (mag [i+62] ) ; /* Last two FFT bins */ 

suma[i+2] = aten* suma [i + 2] + temp; 
wtave = suma [i+2] /denom; 

Unext[i+2] = Uold[i+2] + (temp*temp -2 .0* temp* wtave +wtave * wtave) ; 
ssn = ssn + mag [i+62] *mag [i+62] ; 

> 



/*********** IF past the warmup period, compute the average noise ******** 
************ variance to this time and compute the activity band ******** 

************************* VARIANCE FOR THI S CYCLE ***********************11 

avevar =0; 

i f (k>=warmup) 

{ 

for (i=0; i<4; i++) /* Calculate average outer bin variance 11 

avevar = avevar + 0.25* (k+1) * Unext [i] / (denom*k) ; 

thresh =(klug * avevar); /* Set the threshold */ 

isactive=0; 
ss = O; 
sum ?= O; 

for (i=4; i<60; i++) 

{ 

sss= sss + mag[i] *mag[i] ; /* Calculate activity bin energy */ 

temp = loglO (mag [i] ) ; /*for snr and find var of log mag*/ 

ss = ss + temp * temp; 
sum = sum + temp; 

> 

varib = idf* (ss-sum*sum/ (df+1 .0) ) ; 
if(nar = 1) 

printf ("\ncycle%d varib=%f avevar=%f thresh=%f",k, varib, avevar, threst 

/*************** COMPARE ACTIVITY BAND VARIANCE TO THE THRESHOLD ************ 

if (varib > thresh) 

countO=l ; 
ns = ns + 1; 
if(two=l) 

{ 
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if (countO+countl+count2>=2) /* If threshold is exceeded and */ 

{ /* 2 of last 3 cycles also did so */ 

isactive=l; /* declare activity and write raw */ 

iwrite (x, 128, &irec) ; /* waveform data to output file */ 

> 

> 



else 

{ 

isactive=l; 

print f ("%d active\n" , k) ; 
iwrite (x, 128, &irec) ; 

> 



> 

else 

count0=0; 
nn = nn + 1; 
if(countl + count2 — - 2) 
{ 

isactive=l .0; 
iwrite (x, 128, &irec) ; 

> 

> 

> 

actgraf2 (fp, isactive) ; 

count 2=countl ; 
countl=countO ; 

for (i=0; i<4; i++) 

Uold[i] ' = aten*Unext [i] ; 

k++; 



/* If threshold is not exceeded */ 
/* declare activity only in the */ 
/* event that 2 of the last 3 did*/ 



/* Dump activity status to file */ 
/* for printing later */ 

/* Update counters for 2 out of */ 
/* 3 criterion test */ 

/* Attenuate the numerator of */ 
/* the moving, weighted variance*/ 

/* Count the number of cycles */ 



/******************** CALCULATE SNR EOR THIS CYCLE ********************** 
******** SNR IS THE AVERAGE OVER EACH SIX DECLARATIONS OF ACTIVITY ******/ 

ssn = ssn/4; 
sss = sss/56-ssn; 



if (isactive==l) 

{ 

if (nar— 1) 

printf ("\n**** ACTIVITY ****"); 
if (nisact%6=0) 
avesnr = 0; 
if (sss>0) 

snr = 10.0*logl0 (sss/ssn) ; 
avesnr = avesnr + snr/6; 

> 

if (nisact%6=5 && nar==l) 

printf ("ave snr= %f" , avesnr); 



/* Calculate the average snr */ 
/* over 6 declarations of */ 
/* activity and print out */ 
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nisact++; 

> 

/* Go back to while statement, */ 

} /* go through another cycle */ 

/****************** WHEN WE'RE ALL DONE, PRINT THE NUMBER **************** 
******************* OF CYCLES AND COMFIRM END INDICATOR **************** 

******************* from FILE READING SUBROUTINE ***************/ 



pr int f ( ' 'k=%d end=%d\n ",k, end) ; 
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